Decomposition of the Total Electromagnetic Momentum in a Linear Dielectric into 

Field and Matter Components 
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The long-standing resolution of the Abraham-Minkowski electromagnetic momentum controversy 
is predicated on a decomposition of the total momentum of a closed continuum electrodynamic sys- 
tem into separate field and matter components. Using a microscopic model of a simple linear 
dielectric, we derive Lagrangian equations of motion for the electric dipoles and show that the di- 
electric can be treated as a collection of stationary simple harmonic oscillators that are driven by 
the electric field and produce a polarization field in response. The macroscopic energy and momen- 
tum are defined in terms of the electric, magnetic, and polarization fields that travel through the 
dielectric together as a pulse of electromagnetic radiation. We conclude that both the macroscopic 
energy and the macroscopic momentum are entirely electromagnetic in nature for a simple linear 
dielectric in the absence of significant reflections. 



I. INTRODUCTION 

In a dielectric environment, the momentum of the elec- 
tromagnetic field becomes a contested issue. The matter 
was first raised over a century ago when Minkowski [l|, 
in 1908, and Abraham Q, in 1909, produced conflicting 
expressions for the energy-momentum tensor of an elec- 
tromagnetic field in a linear optical medium. While a 
number of experiments seemed to validate the Abraham 
momentum density 



g. 4 = - (E x H) 



(1) 



other experiments suggested that the Minkowski momen- 
tum density 



g M = - (D x B) 

c 



(2) 



provides the correct description of the momentum of light 
in a linear medium The inability to decide between 
the two energy-momentum tensors and their underlying 
momentum densities became known as the Abraham- 
Minkowski controversy. The resolution of the Abraham- 
Minkowski dilemma was provided by Pcnficld and Haus 
[i| , based on earlier work by M0ller y| , who showed that 
the issue is undecidable because neither the Abraham 
momentum nor the Minkowski momentum is the total 
momentum. Instead, each corresponds to an open de- 
scription of the system in which the momentum is allo- 
cated differently between the field and matter 

Following the work of Penfield and Haus [||, Gordon 
@, in 1973, assumed the Abraham form for the field mo- 
mentum subsystem and derived the material contribu- 
tion to the momentum in a dilute vapor by integrating 
the Lorentz force on each dipole. Gordon derived the 
total momentum density 



for a vapor that is sufficiently dilute that reflections can 
be neglected. Three years later, Peierls 10] followed a 
similar procedure to obtain a total momentum density 



gp 



1 



2c 



(E x B) 



(4) 



that is the arithmetic mean of the Abraham and 
Minkowski forms. The Gordon and Peierls momentum 
densities are equivalent to lowest order in 8n = n — 1. 

The momentum of the field in a linear material con- 
tinues to be a subject of active research. Mansuripur 
flTj [ treats the kinematic momentum of partially reflect- 
ing matter together with the forward traveling light mo- 
mentum, Kemp |12j adopts a subsystem approach for the 
field and matter, and Barnett [13[ shows that the total 
momentum can be regarded, equally, as a combination of 
the Abraham and kinetic momenta or as the sum of the 
Minkowski and canonical momenta. Obukhov [14] stud- 
ied the effect of moving media on the electromagnetic 
energy and momentum. The literature on the subject 
is extensive and the reader is referred to reviews of the 
prior work by Pfeifer, et al. [§[ , Barnett and Loudon [ll| , 
Milonni and Boyd [3, and Kemp [l7j and the resource 
letter by Griffiths [HJ. 

The total momentum of the electromagnetic field and 
dielectric can be uniquely defined only in the context of 
a thermodynamically closed system with complete equa- 
tions of motion Q. In recent work global conserva- 
tion principles were applied to a closed system consisting 
of a quasimonochromatic pulse of electromagnetic radi- 
ation normally incident on a stationary simple dielectric 
with refractive index n through a gradient-index antire- 
flection coating. It was found that the total momentum 



G 



total 



(nE x B) dv 



(5) 



1 



g G = - (nE x B) 



(3) 



is given by Gordon's [9j momentum density, Eq. ([3]), inte- 
grated over all space a. However, Milonni and Boyd [l6| 
revisited Gordon's derivation of the total momentum as 



2 



the sum of the Abraham momentum for the field and 
a material momentum arising from the Lorentz force on 
dipoles and derived the Minkowski momentum, which is 
not conserved, for the total momentum. Conservation of 
the total momentum in a closed system is decisive and the 
Gordon momentum, Eq. ([5]) , is the proven total momen- 
tum for th e lig ht pulse plus stationary dielectric closed 
system [H, Hflf . 

In this article, we investigate the decomposition of 
the total electromagnetic momentum in a linear dielec- 
tric into separate field and matter components. At the 
level of continuum electrodynamics, we take the mate- 
rial momentum to be the difference of the total momen- 
tum, which is conserved, and a field momentum that is to 
be determined. The temporal derivative of the material 
momentum density is the macroscopic force density that 
can be expressed in terms of the fields and a macroscopic 
material property. Assuming that a decomposition of 
the total momentum into field and material components 
exists, we find that each choice of the field momentum 
results in a specific relation between the refractive in- 
dex and the macroscopic material property. Adopting a 
microscopic model of the material as being composed of 
discrete dipoles, we derive Lagrange equations of motion 
and show that the dielectric can be treated as a collec- 
tion of stationary simple harmonic oscillators driven by 
the electric field. The polarization field that is generated 
by the oscillators contributes to both the total energy 
and the total momentum through the refractive index 
n.. We conclude that, in the absence of reflections, the 
total energy and the total momentum are electromag- 
netic in nature and do not separate into field and matter 
components. 



calculations to be performed using the formalism of clas- 
sical continuum electrodynamics. 

We start with a proven fact from the conservation of 
total momentum in a thermodynamically closed system 
consisting of a quasimonochromatic pulse of light in the 
plane-wave limit normally incident on a stationary sim- 
ple dielectric with a gradient-index antireflection coating 
[T9L [20j : The total electromagnetic momentum for the 
case of a negligibly reflecting simple dielectric is given by 
Gordon's formula, Eq. ([5]). The total momentum for the 
system is the sum of its field and material momentum 
components d i, H, i, GJ, 53 ■ Since the 1960s, the reso- 
lution of the Abraham-Minkowski controversy has been 
that any decomposition of the total momentum into a 
field component and a matter component is valid as long 
as their sum is the total momentum [3, 8]. 

The Abraham momentum is a common, but not exclu- 
sive, choice for the momentum of the field in a dielectric 
[H S EH IH| ■ Rather than make this an ansatz, we start 
with a more general form of the field momentum in a 
dielectric 



G field = f ^(ExB)dv, 

J a c 



(6) 



with £(n) an unspecified function of n, that includes 
most of the proposed forms of electromagnetic momen- 
tum, (e.g., the Abraham, Minkowski, and Pierels mo- 
mentums). We note that £ is specifically not a function 
of the electric and magnetic fields. The momentum of 
the material G mat is the difference of the total and field 
momentums, 

G mat = G total — Gf ie id = / (E x B)dv. (7) 



II. FIELD AND MATTER COMPONENTS OF 
THE TOTAL MOMENTUM 

A physical theory is an abstract mathematical model 
of some limited aspect of the real world [2l| , in this case, 
light traveling at speed c/n in a simple linear dielectric. 
Our model dielectric consists of a rectangular prism of 
space in which light travels at speed c/n, where n is a real 
constant. No such idealization of a dielectric exists in the 
real world. Instead, the usual conditions for results of the 
model to correspond to results in the real world apply: A 
simple dielectric is linear, isotropic, homogeneous, trans- 
parent, and dispersion-negligible in a regime in which 
electrostrictive effects and magnetostrictive effects can 
be neglected [22j ■ We assume that a quasimonochromatic 
optical pulse impinges on the dielectric from the vacuum 
at normal incidence in the plane- wave limit and that the 
dielectric is covered with a thin gradient-index coating in 
order to make reflections and surface radiation pressure 
negligible. The material is then stationary in the lab- 
oratory frame of reference and it should be noted that 
the usual Maxwell-Heaviside equations are valid only for 
stationary media. These limits allow our macroscopic 



The material momentum density 

n — ((n) 



(E x B) 



(8) 



leads directly to a macroscopic force density 

_ dg mat . ( dE SB 

f ma t = ~, = (n - C{n j) I ^r— x B + E x 



dt 



d(ct) 



d(ct) 



We define a macroscopic material property 
£ e (n) = n - C(n) 



(9) 



(10) 



that relates the electric and magnetic fields to a macro- 
scopic force density 



, 3E „ <9B 

fmat = £e ( 777 TT X B + E X 



d(ct) 



d(ct) 



(11) 



on the spatially distributed particles of matter. Then 
the decomposition of the total momentum into field and 
matter components, 



G 



field 



n - £ e 



(E x B)dv 



(12) 



3 



G mat = / ^(ExB)*, 

Jo C 



(13) 



is uniquely denned by macroscopic material properties n 
and £ e - 



III. MICROSCOPIC TREATMENT OF 
MATERIAL MOMENTUM 

The total momentum of a closed system with complete 
equations of motion is unique by virtue of its conserva- 
tion. We now know that the decomposition of the total 
momentum into field and material components is not ar- 
bitrary, as had previously been argued H43 [HHH] ■ 
We simply need to determine the nature of the material 
momentum and what those properties are. 

The most rigorous approach to light propagation in a 
material is to use a microscopic model in which the mi- 
croscopic electromagnetic field interacts with individual 
electric dipoles in the vacuum. If the distance between 
dipoles in the material is much smaller than the wave- 
length, then the theory can be considerably simplified 
by applying an averaging scheme to a volume contain- 
ing many atoms per cubic radiation wavelength. Gordon 
0, for example, spatially averages the Lorentz force on 
a single electric dipole with dipole moment p = ae 



fdipoie = (P • V)e + x b 

a(ct) 

to create a macroscopic Lorentz force density 



(14) 



N(f t 



dipole/ 



Na 



E x 



<9B 



<9E 



d(ct) d(ct) 



x B 



(15) 

Here, e is the microscopic electric field, b is the micro- 
scopic magnetic field, a is the linear polarizability, N is 
the number density, E = (e) is the macroscopic electric 
field, and B = (b) is the macroscopic magnetic field. 
Gordon Q employs the low-density approximation 



Vl + Na w 1 + Na/2 



(16) 



such that 



(17) 

Other authors have derived variations of this basic rela- 
tion. Peierls [13 used the exact relation 



n = 1 + Na 



(18) 



to obtain 



, n 2 - 1 / dB <9E 



Saldanha 0| and Bradshaw, et al. [H[ find the macro- 
scopic Lorentz force 

(20) 

If the Lorentz force, Eq. ([20]) . is separated into compo- 
nents, 



fza = in — 1) I E x 



dB 



<9E 



d(ct) d(ct) 



x B 



(21) 



(22) 



then f^2 is in the form of the Abraham force and differs 
from the Peierl's result for the Lorentz force by a factor 
of 2. Selecting one of these average Lorentz force den- 
sities to equate with the macroscopic force density f m at, 
Eq. (|lip , will result in a value for the material parameter 
£ e and formulas for the field and material components of 
momentum in a dielectric, Eqs. (p~2]) and (fl~3|) . 

The Gordon form of the force, Eq. (|17|). provides us 
with the sensible result that the material momentum is 
the Abraham momentum. But the Gordon form is an ap- 
proximation. Peierls's exact result, Eq. (fl"9|) . provides us 
with a field momentum that is negative for some values 
of the refractive index. Further, there is a factor of two 
discrepancy in the derivations of Gordon |9J and Peierls 
[lol ] compared to the later work 0, 23] . The Lorentz force 
given by Eq. (|22|) presents us again with the problem of 
a negative field momentum for some values of n. There- 
fore, none of the versions of the Lorentz force that are 
presented here are particularly appealing. In the next 
section, we take a more detailed look at the microscopic 
Lorentz force. 



IV. LAGRANGIAN DYNAMICS OF THE 
DIELECTRIC 

The macroscopic force density, or Lorentz force, is only 
significant where the amplitudes of fields are changing in 
time. The effect of the Lorentz force has been inter- 
preted, in the context of a "top hat" pulse shape, as an 
acceleration of the physical atoms/dipoles on the front 
side of the pulse where the electromagnetic field is in- 
creasing in strength 9]. The atoms in the middle por- 
tion of the pulse travel at a generally constant velocity 
until they are decelerated at the back edge of the field. 
No momentum is left in the atoms after the pulse has 
passed. The process has been described by saying that 
the material momentum travels with the pulse [9J . 

The Lagrangian for a collection of electric dipoles in- 
teracting with a microscopic electromagnetic field is [HI 



(19) 



- -KY 2 r 

2 c 



1 



■MR; 



1 



4 



+ ^ f Q -A(R Q ) 



M 



r a • [(r Q • V)A(R Q ) 



1 

2? 



+R Q • [(r„ • V)A(R Q )] 



(A(x)) 2 -c 2 (V x A(x))' 



d 3 x. 



(23) 



where a labels each dipole, M = m+ + m_ is the total 
mass, n = m+m_/M is the reduced mass, r = r + — r is 
the relative coordinate, R = (m + r + + m_r_)/M is the 
center-of-mass coordinate, and A is the vector potential. 
For each a, the Lagrange equations of motion are 



d ( dL 



dL 



dt KdRiJ dR 



d f dL 
dt V dii 



dL_ 



= 0. 



(24) 



(25) 



Applying Eq. (|M|) to the Lagrangian, Eq. ([2"3"]l and drop- 
ping the less significant terms, we obtain a microscopic 
force equation 



MR = -2(r • V)^ + -r x (V x A) 

c at c 

that is the Lorentz force on a single dipole, Eq. 
The other Lagrange equation of motion 



(26) 



(j,r + kv 



-(R- V)A 



(27) 



that is derived by applying Eq. (|25l) to the Lagrangian, 
Eq. (|23p , and keeping the significant terms is more in- 
teresting because it is evocative of the other microscopic 
model of a dielectric as stationary simple harmonic oscil- 
lators driven by the electric field. 

Momentum is a quantity that represents the movement 
of energy. Typically, this is mass energy and it is nat- 
ural to think of the atoms that are accelerated by the 
Lorentz dipole force as having momentum p = mv. In 
this Gordon [9j model, the atoms are free particles that 
are accelerated by the Lorentz force at the leading edge 
of the field and travel at constant velocity until decel- 
erated by the Lorentz force at the trailing edge of the 
pulse. Then, the momentum of the atoms contributes to 
the total momentum. However, that model assumes a 
low-density limit in which there is no scattering of the 
field and assumes that the atoms move without being 
restrained by association with a lattice and without col- 
liding with any other particles in the material. That is 



not the case for a dielectric and such a material has no 
definite index of refraction. Instead, there is substan- 
tial impediment to the motion of atoms in a dielectric 
that, after all, is a material that can be modeled as a 
continuous medium. By setting R = in Eq. fl27[) we re- 
cover the usual and familiar stationary simple harmonic 
oscillator model of a dielectric. The oscillators with fre- 
quency loq = \J k/h, driven by the quasimonochromatic 
field with center frequency w p , create a polarization field 



E 



(28) 



that contributes to both thetotal energy and the total 
momentum through the refractive index n. The fields, 
electric, magnetic, and polarization, travel together as a 
pulse of electromagnetic radiation and are responsible for 
the electromagnetic momentum. Although there will be 
some small and transient amounts of energy and momen- 
tum associated with the action of the Lorentz force on 
massive particles, the energy and momentum are almost 
entirely electromagnetic in nature. Then the momen- 
tum density of the matter, Eq. (J5J, is negligible, as is 
the associated macroscopic force density, Eq. (j9|). Sub- 
stituting £ e = into Eqs. (IT21 and (llUl) . we have and 
Gfield = Gtotai and G mat = 0. There is no apprecia- 
ble matter component of either energy or momentum for 
a simple linear dielectric medium that is at rest in the 
laboratory frame of reference. 

The situation changes somewhat in the absence of a 
gradient-index antireflection coating. Momentum bal- 
ance requires, at the cost of some amount of rigor in 
the derivation, that we impute a momentum to the ma- 
terial that is twice the momentum of the reflected field 
imparted by the surface force of reflection. But, that is 
the only significant material momentum, there is no mo- 
mentum of the Gordon type due to motion of the atoms 
inside the dielectric. 



V. CONCLUSIONS 

The total momentum is uniquely conserved in a ther- 
modynamically closed system with complete equations 
of motion. Once the total momentum is known, with 
certainty, from conservation principles, then the momen- 
tum relations can be derived. Using a microscopic model 
of the dielectric medium, we showed that the propagat- 
ing polarization field, not the movement of matter, is a 
component of the momentum of a system consisting of 
a stationary dielectric illuminated by a pulse of quasi- 
monochromatic radiation through a gradient-index an- 
tireflection coating. 
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